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Q\ ■ Abstract 



The Jacobi evolution method has been widely used in the QCD analysis 
of structure function data. However a recent paper claims that there are 
serious problems with its convergence and stability. Here we briefly review 
the evidence for the adequate convergence of the method; and show that there 
are errors in the above paper which undermine its conclusions. 
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In analysing data on nucleon structure functions, it is convenient to have as fast and sim- 
ple a method of implementing their QCD evolution as possible. For this reason, the Jacobi 
evolution method |L]-|3| has been adopted in many recent analyses^, including applications to 
polarized [|| and non-forward or off-diagonal |JJ parton distributions. However, Ghosh and 
Raha have made the startling claim that it is unstable due to convergence problems. In 
this paper, we discuss the convergence properties of the method in general, and the problem 
raised by Ghosh and Raha in particular. 

We start by very briefly summarizing the method |T|-j3|. The basic idea is to expand a 
given structure function A(x, Q 2 ) at fixed Q 2 in the form 



A(x, Q 2 ) =xf>(l- x) a £ a k (Q 2 )ef(x) , (1) 



k=0 



where a, (3 > — 1 are real numbers and 0^ (x) are a complete set of Jacobi polynomials of 
order k — 0, 1, 2, They are defined to satisfy^ 

dxx?(l - x) a &f(x) ef(x) = 5 kl , (2) 







A SPIRES search reveals more than 30 citations since 1995. 



2 We follow the standard notation and conventions of J2|, where further mathematical details may 
be found. 
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so that the "Jacobi moments" 

a k (Q 2 ) = f 1 dxA(x,Q 2 )Qf (x) . (3) 
Jo 

In practice, the truncated expansion 

N 

A(x, Q 2 ) = A N (x, Q 2 ) = x?(l - x) a £ a h (Q 2 )ef(x) (4) 

k=0 

is used, where N is hopefully large enough for this to be a good approximation. The structure 
function can then be evolved in Q 2 by using the appropriate QCD evolution equations for 
the "Jacobi moments." For non-singlet structure functions these are of the form 

a k (Q 2 ) = j:E^(Q 2 ,Q 2 ,A)a 1 (Q 2 )i (5) 

3=0 

where the £jf(Q 2 , Qg, A) are known coefficients and Qq is some conveniently chosen ref- 
erence Q 2 . Alternatively and equivalently, the Jacobi moments dfc(Q 2 ) can be expressed 
in terms of the more familiar Cornwall- Norton moments at the reference Qq. For singlet 
structure functions, the evolution equations are modified by the addition of appropriate in- 
homogeneous terms dependent on the moments of the gluon distributions at the reference 

Ql 



II. CONVERGENCE OF THE JACOBI EXPANSION 

The conditions for uniform convergence |3||| of the expansion (|T|) in the range < x < 1 
are the existence of the integrals 

f 1 dxx^l -x) a A(x,Q 2 ) [ dxx^l-x) 

Jo Jo 

for < z < 1. For non-singlet structure functions, assumed to be analytic and bounded, 
these are satisfied for the whole range of values a, (3 > —1 for which Jacobi polynomials 
exist. For singlet structure functions, assumed to behave like x~ 7 as x — > with 7 > 0, 
they are satisfied for a > — 1, (3 > — 1 + 27. The key question is not whether a convergent 
series exists, but whether it converges sufficiently rapidly for the truncated expansion 
to be a good approximation for suitably chosen values of a and (3. These values can reflect 
the observed behaviour at small and large x, or simply be varied in order to improve the 
convergence. 

The accuracy of the truncated expansion (Q) can be checked in several different ways. By 
construction (|4]|J) automatically guarantees the correct evolution of the first iV Cornwall- 
Norton moments, but other QCD evolution properties will only be correctedly reproduced 
if the series (f|)has converged. For example, integrals of the form 

f dxA(x,Q 2 ) I x (7) 
J 



A(x,Q 2 )-A(z } Q 2 ) 



x — z 



(6) 
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are sensitive to the \ow-x region, and correspond to the well-known Gross-Llewellyn-Smith 
and Gottfried sum rule integrals for A = xF 3 and A = F 2 ep — F 2 era respectively. These are 
independent of Q 2 in LO QCD, but this is only guaranteed by (|j) if the series has converged. 
Barker, Langensiepen and Shaw checked this property in their LO fits and found that 
with N — 10 it was satisfied over the range 2 < Q 2 < 200 GeV 2 to within an accuracy of 
10 -3 . They also checked that the QCD evolution of the higher Cornwall-Norton moments 
Mi(Q 2 ) for i = 15, 20 was accurately reproduced?] and concluded that truncating at N = 10 
was sufficient for all practical purposes over the kinematic range considered. 

Other methods are to check the variation of the predictions of (|j, |5|) with increasing 
iV; and to compare them with the results of "exact" numerical evolution of the DGLAP 
equations for the same input. Both methods were exploited for non-singlet strucure functions 
by Chyla and Rames J7|, who compared the results of "exact" QCD evolution with the 
results obtained by trucating the Jacobi series at N values varying up to 20. Specifically, 
they assumed the input form 

A(x,Q 2 ) =4.05x°- 404 (l-x) 3 - 73 (Qo = 90) (8) 

and considered its evolution over the Q 2 range 1 < Q 2 < 1000 GeV 2 . They found a rapidly 
convergent oscillatory approach to the "exact solution" using the same values a = 3, (3 = 1/2 
as [Q]; and that small variations in these values had little effect on the convergence. For 
example, detailed results are presented for the reasonable choice Klo — 0.1, retaining terms 
up to N = 9 . Even at x = 0.7 where the structure function is very small, it is reproduced 
with an accuracy of < 1%, corresponding to an absolute error of < 0.0005 This is very small 
on the scale of variation of the structure function itself, which has a maximum value of order 
unity. 

Chyla and Rames |7j] also made similar studies for the other evolution methods which 
had been proposed and concluded that the Jacobi method was "by far the simplest, fastest 
and simultaneously very accurate way of analysing the Q 2 evolution of structure functions. " 

Finally, we note that for singlet structure functions, Krivokhizhin et al. H have presented 
results for the errors obtained in reproducing a "typical structure function" 

A(x) = 2.67x a25 (l - x) 3 + 0.48(1 - x) 8 

at fixed Q 2 using the truncated expansion (U) for N values up to 17. For the most favourable 
values a = 3, j3 = —0.85 they found an average error < 10~ 4 for N = 7; and for N = 8, the 
error at any x was always less than 10 -4 . For other a, f3 values the convergence is less rapid 
but still more than adequate for practical purposes. 

III. A POSSIBLE PROBLEM? 

Very recently, Ghosh and Raha have again investigated the convergence of the Ja- 
cobi method by assuming specific analytic forms for non-singlet structure functions at the 



3 In LO, the logarithm of the Cornwall-Norton moments depends linearly on tin [In Q 2 /A 2 ] . This 
linear dependence was reproduced with a fitted slope within one per cent of the QCD prediction. 
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reference Qq. These are approximated by the expansion (||) using (|3|) to evaluate the Jacobi 
moments at Qq and @ to evolve them to other Q 2 . 

The basic problem is already manifest at the reference Q 2 = Qq. Plots are presented 
showing that as N is increased from zero, their computed value of the truncated series 
An(x, Qq) initially oscillates and then becomes approximately stable and equal to the exact 
A(x, Qq) for a wide range of iV values. This is similar to the behaviour found by other 
authors. However when N increases further, beyond a value which is typically of order 20 - 
25, the results become wildly unstable, presumably reflecting large contributions from high 
order polynomials in the Jacobi series ([!]). They therefore conclude that "the convergence 
breaks down completely for large N. " 

This reult is surprising, since the input structure functions are smoothly varying and the 
series must eventually converge uniformly, as noted above. However the precise evaluation of 
both Jacobi polynomials and of the integrals (§) becomes non-trivial as the order increases; 
and no tests are reported with a view to verifying the accuracy of the methods used for 
the highest orders considered. The question therefore arises whether the reported problem 
is really endemic to the Jacobi method; or whether it could arise from problems with the 
procedures used, which are not specified in any detail. Fortunately, this can be easily settled 
using the first set of input forms considered, which are of the generic typeQ 

A(x, Ql) = x" 2 {\ - x 2 f f>(l- x 2 ) k . (9) 

fc=0 

This can obviously be rewritten in the form 

2p+4 2p+4 

A(x,Q 2 Q ) = x^il - xf £ c, ^ . = x^{l - xf £ aj {Ql)Qf /2 {x) , (10) 

j=i j=l 

where we have assumed the usual values a = 3 j3 = 1/2 for the Jacobi indices. In other 
words, ak(Ql) = for k > N max = 2p+4 and the Jacobi expansion, implemented in the most 
straightforward way, terminates exactly after a small number of terms. The problematic 
contributions found by Ghosh and SitaiQ at large iV are completely excluded in this case. 
They are not endemic in the Jacobi method, and must arise from some problem with their 
procedure, numerical or otherwise. 



4 Cf. eqns.(ll) of for A = xu v (x), xd v (x) 

5 Another problem raised in their paper is that even if one restricts oneself to N values in the 
range where the results are stable, their results are strongly dependent on whether they start their 
evolution from Q 2 = 3.5 or Q 2 = 5 GeV 2 . This stems from two trivial errors in the equations used 
to describe the input at Q 2 = 3.5. For full details see the earlier and longer version of this paper, 
which is still accessible on lanl hep-ph/9901253| vl. 
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IV. SUMMARY 

Previous work has established that the Jacobi expansion coverges very rapidly for rea- 
sonably chosen parameters a, f3. In the light of the previous subsection, we do not believe 
that Ghosh and Raha have provided any reason to modify this conclusion. 
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